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Abstract 



We compute non-perturbative contributions to B — > X s l + l~ that are not 
explicitly suppressed by powers of the 6-quark mass. They are proportional to 
(B\ha-Gh\B) and arise from an interference between the free-quark amplitude 
and higher order terms in the matrix element of the four-quark operator 
s Q! 7 At (l — 7s)c Q c^7 M (l — 75)6/3. This correction is found to be small over 
most of the dalitz plot except near the charm threshold. Unfortunately, the 
perturbative computation we have performed is invalid near charm threshold 
and we do not except to see the structure found at lowest order reproduced 
in the data. We conclude that these non-perturbative contributions do not 
significantly modify the previous analysis of B — ► X s . 
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Rare decays of B mesons such as B — > X s 7 and B — > X s provide a probe of physics 
at distance scales much smaller than the compton wavelength of the b-quark. The rates 
for these processes vanish at tree-level in the standard model of electroweak interactions 
but proceed at one-loop and beyond. Their GIM "suppressed" amplitudes are dominated 
by the loop contributions from the top quark. As such, B — > Xgj and B — >• X s are 
ideal candidates for observing physics beyond that of the standard model, e.g. an extended 
Higgs sector. However, in order to use them to identify or even claim a hint of new physics 
one must have a robust prediction for decay rates and/or other observables in the standard 
model alone. Both the inclusive B — ► X s 7 (with cuts on the photon energy E^) and some 
exclusive radiative decay modes such as B —>■ K*^ have been observed. As yet neither 
inclusive nor exclusive decays B — > X s have been observed, however, the CLEO Q 
upper limit of 1.6 x 10~ 5 on the B° — ► K*°e + e~ branching ratio is within a factor of ~ 3 
of a standard model estimate. As theoretical predictions for exclusive decay widths are 
notoriously difficult to obtain it is the inclusive decays that one is presently looking toward 
to make a meaningful comparison between theory and experiment. 

The heavy quark limit of QCD leads to great simplification in the theoretical analysis 
of inclusive decays of B mesons [0]. In particular, the inclusive decays of B mesons can 
be described using an operator-product expansion since the characteristic momentum of 
the process ~ m& is large compared to the scale of strong interactions ~ Aqcd- A sys- 
tematic expansion of decay observables, such as the rate and differential spectra can be 
performed in powers of Aqcd/^6 and the strong coupling constant a s (fj). There are three 
non-perturbative matrix elements, A, Ai and A2, that appear in expressions for such 
observables up to order I /ml in the expansion. 

The analysis of the inclusive decays B — > X s ^y and B —> X s l + l~ [|] to O (1/mf) is 
slightly different to the analogous computations performed for inclusive semileptonic decays. 
The difference lies in the fact that the semileptonic decays proceed at tree-level from charged 
current interactions while B — > X s ^ and B — > X s l + l~ first occur at one-loop. There are 
contributions to B — > X s ^ and B — > X s from closed charm quark loops in addition 
to local operators in the effective theory resulting from integrating out the top quark and 
W ± , Z°-bosons (we do not discuss up-quark contributions). It was recently suggested that 
these loop contributions give power corrections to the rate for B — > X s ^ suppressed by 
the mass of the charm quark |J as opposed to the mass of the bottom quark @,f|. Explicit 
computation || showed that these new contributions are proportional to the chromomagnetic 
matrix element (B\ha ■ Gh\B). Subsequently, it was shown that such contributions 
to B — > X s 7 do not dominate the I /ml power corrections found in In addition 

it was pointed out [7|,j8| that there is an infinite tower of operators whose contributions 
are suppressed only by powers of m\,KQCj)/m 2 c but which appear with numerically small 
coefficients. 

The Feynman diagrams highlighted by Voloshin yield a term that vanishes for on- 
shell photons in addition to the term modifying B — > A s 7. In this work we determine 
the contribution of these diagrams to the inclusive decay B — > X s where I denotes 
an electron or a muon. We do not re-examine the O (1/mf) contributions [|],|9] as we have 
nothing new to add. Further, we do not include the non-renormalization group pQCD 
corrections to the decay width ]nj. Only the leading free-quark amplitude and the loop 
correction yielding this new non-perturbative contribution are considered. 
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For momentum scales much below the mass of the W ± gauge boson it is appropriate to 
describe these rare decays by an effective hamiltonian 



AG 



V 2 i 



(1) 



where the are Wilson coefficients of the operators Oi(fj) renormalized at a scale \i 

TT|— p~3| . We have assumed that VuV} t = — HcKc i n writing T-C e s this way and we assume 



that the leptons and strange quark are massless. In order to fix notation we list the operators 
(we use the same notation as |TT| ) 



o 2 
o 3 
o A 
o 5 
o 6 
o 7 

o 9 

Oxo 



SL,al^b Lja Ul,P7ij.Ul,I3 + 
SL,on^b Lj p UL,j3lnU L ,a + 

SL,aYb L) p Ur^Ur^ + 
g 

m b s L:Q a^ u b R 
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■^R,pliJ,bR,a 
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e 2 
I6tt 2 

e 2 

I6tt 2 
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mgiL^KpbRs a 



a 

{3 ^fiu 
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Working at leading order |TT[ we truncate the operator basis to Ox, 2 , 7 , O s and Og. 
The next to leading order results can be found in |T2| , |T3"| and a complete summary including 
the non-perturbative and pQCD corrections can be found in |J. 

The rate and other "short-distance" observables for the decay arise from an operator- 
product expansion of 



T = Im 



d*x e~ iqx {B\ T H eS {x)nl s (0) \B) 



(3) 



integrated over the appropriate leptonic phase space. The free quark decay result emerges 
as the leading term in the OPE of T with higher order corrections to the free quark result 
arising systematically £|,[J as a power series in and 1/m". The component of T 

previously omitted from B — > X s ^ arises in 



T 



c 2 2 (x) c* 7 O\(0) + h.c. 



(4) 



where the charm loop closes off via a soft gluon and a photon, as shown in figure 1. The 
diagrams that modify B — > X s l + l~ are shown in figure 2 (we have not shown the hermitian 
conjugate). 
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The amplitude resulting from closing the charm quark legs of 2 with the emission of a 
gluon of momentum k with polarization 7], and a photon of momentum p with polarization 
A is given by (to leading order in k) 



A = ^ sTV(l - T5)& 



h \v^p a k h ~h\p- ke a ^p a + I 2 p 2 e Xb ^k b - I 2 Px e ab ^p a k b 



where the functions 7^2 are given by 

Ii2= f 1 dx r 1 '* d 2x y 1 / 6 n 

io •/ [m 2 . — x(l — x)p 2 — 2xyp ■ k] 

As k 2 ~ Aq CD <C m 2 . we neglected it in both the numerator of A and denominator of Ii j2 - 
However the p-k term appearing in the denominator of Ii j2 is of order ~ E 1 Aqcd ~ m 2 and 
dominates over the other terms near the charm pair threshold (p 2 ~ 4m 2 ) region. Before 
proceeding it is interesting to explore the issues raised in |7]||. We expand the integrals Ix t2 
in terms of p-k and evalute the coefficients of terms higher order in Aqcd- After performing 
the angular integration over phase space one finds that 

00 

<I 1 (s,m 2 c )> = a nK{s) (J5 7 Aqcd)" 

n=even 

1 / \ f 1 , x n+1 (l-x) n+2 

K(s) = dx— +1 , (7) 

Jo [m A c — sxyl — X)) 

where E 1 is the energy of the lepton pair in the rest frame of the B-meson. The first few 
constants a n are ao = \ , = I , 04 = j|. In the case of on-shell photons, s = p 2 = 
and E 1 = one finds that the contribution from terms higher order in Aqcd converges 

rapidly. When the photon is off its mass-shell as is the case for B —>■ X s and the 
invariant mass of the lepton pair are far from the charm threshold region the series is again 
convergent (requiring that m^ — \fs ^> Aqcd)- In contrast, when the invariant mass of the 
lepton pair is near the charm threshold region the series is not convergent. This is to be 
expected since we know that the charmonium resonances in this region do not appear at any 
finite order in perturbation theory. In fact, this problem is not isolated to our computation 
but is a problem for the entire parton level computation of the differential rate. One expects 
that suitable smoothed quantities over this region of invariant mass will be reproduced by 
the same smoothed quantities in the parton model by duality. For the rest of our discussion 
we will drop the p ■ k terms so that the integrals 7^2 become identical 1\ = I2 = I where 

H*,ml)= f ^ , 2 V{1 ~ V) \ ] ■ (8) 
Jo [mi - y(l - y)s] 

The leading contribution to the differential decay rate comes from the free quark decay 
and at leading order in strong interactions [I1J is 



jp(0) 2 f<2 5 

' r " b \V bt V s \\ 2 (1-s) 2 (l + 2s)\c 8 (m b )+g(r,s)(c 2 (m b ) + 3c 1 (m b ))\ 2 



ds An 2 192tt j 



(5) 



+ (l + 2s)|c 9 (m b )| 2 + 4(l + 2/s)|c 7 (m b )| 2 
+12Re [c 7 {m b )*{c 8 {m b ) + g{r, s){c 2 {m b ) + 3 Cl (m 6 )))]] , (9) 



4 



where r = m c /mb is the ratio of quark masses and s = p 2 /m\ is the invariant mass of the 
lepton pair. The function g(z,s) arises from closing the charm loop from operators (9i j2 
forming a virtual photon, and is given by (at fi = m^) |Tl|-p^] 



-—Inz H 1 — v 

9 27 r 



-l(2 + y)J\l-y\ 



6(1 -y) [in 



1 + 



+ m\ +6{y - 1)2 tan- 1 



■ (10) 



4z 2 



with y 

Computation of the non-perturbative term is straightforward but tedious and we will 
not detail it here. The additional contribution to the differential rate is found to be 



ds 



a 



Glml 



\V bt V ] st \ 2 



32A 5 



3ml 



4tt 2 192tt 3 
Re [( (1/s + 6 — s) Ci^rrib 
c 2 (m 6 ) I(s,r) ] 



+ (2 + s) (c 8 (m 6 ) + (c 2 (m b ) + 3 Ci(m fe )) g-(r, s) ) )* 



The nonperturbative matrix element A 2 that appears in this contribution is found from 
the B* — B mass difference to be A 2 ~ O.lGeV 2 . Comparing the size of this new non- 
perturbative contribution to the differential decay rate, as shown in figure 3, one see that it 
is small except near the charm threshold. Integrating the expression for the rates between 



s = ml/2 and s 



m, 



we find the rate to be decreased by ~ 1.5% due to this new contribution 



(i.e. smaller than the contribution from the I /ml operators, ~ 3%). Unfortunately near 
threshold is precisely the region where the above computation fails due to the presense of 
incalculable long-distance physics associated with the charmonium resonances. In the region 
either far above or far below the charm threshold the correction considered here can be safely 
neglected. 

In conclusion, we have determined a previously omitted non-perturbative contribution 
to B — >• X s , that is not explicitly suppressed by 1/m^. This has the same origin as the 
1/ml contribution to B — > recently observed by Voloshin |J. We find that in the regions 
of lepton invariant mass where we can reliably compute such contributions they are small. 
The computed amplitude has interesting threshold behavior near s = 4ml but unfortunately 
large strong interactions corrections that we cannot determine will significantly modify the 
result in this region. 



ACKNOWLEDGMENTS 

This work is partially supported by the US Department of Energy. 



5 



o 2 o 7 



FIG. 1. A contribution to B — ► X s j from the time-ordered product T[ Tt e s(x) 7i cS (0) ] pro- 
portional to A2 and naively suppressed by 1/mji. 
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FIG. 2. Contributions to B — > X s l + l from the time-ordered product T[ 7i c s(x) H 
proportional to A2. 




FIG. 3. The differential decay rate for inclusive B — > X s l + l~ using quark masses m c = 1.5 GeV 
and mfe = 4.5 GeV normalized to the semileptonic width. The dashed line is the leading order, 
free-quark contribution and the solid line is the sum of the free-quark and the non-perturbative 
contribution computed in this work. Inset is the distribution in the charm threshold region. We 
have used c\ = 0.22, C2 = —1.09, = 0.31, cs = —4.16 and eg = 4.87 at the scale fi = m&. 
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